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Introduction: We present a comprehensive 

visible/near- infrared two-layer radiative transfer 
modeling study using laboratory spectra of variable 
dust thicknesses deposited on substrates with known 
photometric parameters. The masking effects of 
Martian airfall dust deposition on rocks, soils, and 
lander/rover components provides the incentive to 
improve two-layer models [1-3]. It is believed that 
the model presented will facilitate understanding of 
the spectral and compositional properties of both the 
dust layer and substrate material, and allow for better 
compensation for dust deposition. 

Model: We have implemented an adaptation of 
the Hapke model of bidirectional reflectance of a 
two-layer medium ([4] p. 251). This adapted model 
allows the particulate lower layer in the two-layer 
model to be replaced with an arbitrary substrate 
defined only by its Bidirectional Reflectance 
Distribution Function (BRDF). The freedom of 
definition for the substrate material allows for the 
accurate modeling of dust accumulation on non- 
Hapke materials (e.g., rocks with a strong specular 
scattering lobe, silicone rubber RTV used in 
calibration targets). 

This adaptation of the Hapke two- layer model 
consists of the substitution of the substrate’s single 
particle angular scattering function p L (g ), which 
depends only on phase angle and is meaningful only 
in the context of particulate media, with an analogous 
bidirectional scattering function qiiz,el,az',eV), where 
az,el define the incident vector, and az',el' define the 
emission vector. This bidirectional scattering 
function ( q ) is defined as: 
q{az,el,az\eV) = jt • BRDF (az,el,az\el } )/ r s 
where r s is the spherical reflectance. In practice, the 
spherical reflectance is frequently unknown, and 
normalization of q is performed by numerical 
integration of the BRDF. 

The bidirectional scattering function ( q ) 
simplifies to the single particle angular scattering 
function (p ) in the case of a Lambertian substrate. In 
addition the bidirectional scattering function fulfills 
an analogue of the single particle angular scattering 
function’s normalization constraint: 

J J'q(az,e/,az\e/')dQdQ' = (2jt) 2 

2 jt 2 jt 

The Hapke two-layer model also depends on the 


substrate albedo factor (g L ). Although g L can in 
theory be derived from the BRDF 1 , concern over 
error implicit in the calculation of r 0 from r s and over 
the potentially magnified effect of errors in the 
substrate BRDF led us to numerically fit rather than 
analytically solve for g L . 

The specific substrate model used was that 
developed by the Mars Exploration Rover (MER) 
Panoramic Camera team to describe the Panoramic 
Camera Radiometric Calibration Target (RCT) [6]. 
This model consists of the He-Torrance model [5] - a 
physical optics model borrowed from the realm of 
computer science - combined with a Hapke 
backscatter term [4]. A three-parameter Henyey- 
Greenstein function was used to fit the upper layer 
(dust) phase function. 

Data: The Bloomsburg University Goniometer 

(BUG) was used to acquire bidirectional reflectances 
of the MER Pancam RCT materials (silicone rubber 
RTV surfaces with approximately 20%, 40%, and 
60% reflectances in the visible/near-infrared) at four 
wavelengths (480, 600, 750, and 930 nm). 
Measurements also were acquired with variable mean 
thicknesses (0 to 225 pm) of Mars analog JSC-1 dust 
deposited on these substrates using an airfall settling 
technique [7]. 

Procedure: The bidirectional two- layer reflect- 
ance model was fit to the entire data set using a 
Levenberg-Marquardt least- squares minimization 
routine with numerically calculated derivatives. 

The numerical integration of the substrate models 
for purposes of normalizing q L was restricted to 30- 
50° elevation for the emission vector and 20-90° 
elevation for the incidence vector. This prevented 
poor characterization of the RCT BRDF for vectors 
well outside those acquired using the BUG from 
negatively impacting the performance of the two- 
layer model. 

Results/Discussion: The model fit the data with 
a reduced chi-square of 9.1 when the BUG data was 
assumed to possess a relative error of 5%. The 
efficacy of this fit can be visually judged in Figures 1 


Spherical reflectance (r s ) can be calculated directly 
from the BRDF, diffusive reflectance (r 0 ) can be calculated 
from r s ([4] p. 269), the volume single scattering albedo 
(w L ) can be calculated from r 0 ([4] p. 291), and the albedo 
factor (g L ) is defined in terms of w L . 
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Thickness 


0 um 
5 um 
10 um 
24 um 
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132 um 
225 um 


Wavelength 
4 80 nm 
600 nm 
750 nm 
930 nm 
4 80 nm 
600 nm 
750 nm 
930 nm 
4 80 nm 
600 nm 
750 nm 
930 nm 


Value 


All Data Points 


0.16534794 

0.054071909 

0.033830598 

0.042430620 

-0.75249348 

-0.84302023 

-0.86707857 

-0.87971395 


4 80 nm 
600 nm 
750 nm 
930 nm 
4 80 nm 
600 nm 
750 nm 
930 nm 
4 80 nm 
600 nm 
750 nm 
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4 80 nm 
600 nm 
750 nm 
930 nm 


97656175 

99224375 

99538410 

99653797 

99999e-05 

99999e-05 

008792356 

038402812 

31140855 

1573110 

3772483 

59818068 

55086283 

48988512 

43596463 

83078440 

84260554 

79792799 

74477210 

88270952 

91779216 

89677571 

88219746 

49142599 

85253148 

96346270 

96832510 


Table 1 - Parameters fit by the model and their associated values. Where 
parameters depend on substrate material, dust thickness, or filter 
wavelength that information is also provided. PuO, Pul, and Pu2 are 
respectively the forward asymmetry, backward asymmetry, and forward 
fraction parameters for the upper layer three-parameter Henyey- 
Greenstein phase function, tau is the dust opacity, and wu and wl are the 
single scattering albedos of the upper and lower layer, respectively. It 
should be noted that the near unity forward fractions (. Pu2 ) make the large 
backward asymmetry parameters (Pul) nearly meaningless, as the 
backward lobe of the upper layer phase function is nearly nonexistent. 
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Figure 2 - Measured (closed symbols) and modeled (open symbols) spectra 
for several dust thickness for gray (40%) and black (20%) substrates. 



Figure 1 - Scatter plot of measured vs. fit BRDF values 
for all substrates, dust thicknesses, and geometries. Red 
line is 1:1 correlation line. Data points from the white 
(60%) substrate are highlighted green. Data collection 
was hampered for the white substrate due to the small 
size of the available silicone RTV sample. 

and 2, which show a scatter plot of measured vs. 
modeled data, and plots of measured and modeled 
spectra (BRDF * n) for different substrates and dust 
thicknesses respectively. 

The data acquired of the white (60% reflective) 
RCT substrate is marked green in Figure 1. The 
deviation of these measured values from the model is 
consistent with a known data acquisition problem 
relating to the size of the white silicone RTV sample. 
Indeed, when the fit is run without the white substrate 
data the reduced chi-square takes on the improved 
value of 7.6. 

The parameter values fit by the model (see Table 
1) show physically realistic trends. Fitted optical 
depth tau correlates well with measured coating 
thickness, and derived single particle phase functions 
for the dust are nearly isotropic. Derived dust single 
scattering albedos exhibit the expected red slope. 
Figure 2 shows that the model does a good job of 
matching the spectral evolution of the surface as 
more dust is deposited, with the gray and black 
substrates being progressively darkened at blue 
wavelengths, and progressively brightened at red 
wavelengths. 
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